AD-A066  410  MOORE  SCHOOL  OF  ELECTRICAL  EN6INEERIN6  PHILADELPHIA  PA  F/6  12/1 

ROBUST  WIENER  FILTERS  FOR  CORRELATED  SIGNALS' AND  NOISE. tU> 

1960  G  MOUSTAKIDES.  S  A  K ASS AN  AFOSR-77-315* 

AFOSR-TR-80-0426  NL 


UNCLASSIFIED 


ADA  0  8  64 10 


Stfi 


UNCLASSIFIED 


£ 

.  < 

i* 


f  ► 

*, 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (WlMl  Data  F.nli-rvd) 


i 


u~ 

c: 


1.  REPORT  NUMBER  /  2.  GOVT  ACCESSION  NC 

AFQ5R-1R-  8  0-0426  \ fr.fo  T< 


REPORT  DOCUMENTATION  PAGE 

GOVT  ACCESSION  NO 


4.  TITLE  fane/  5u6/ll/eJ  « 

ROBUST  WIENER  FILTERS  FOR  CORRELATED  SIGNALS 
AND  NOISE 


7.  AUTHORfs; 

G.  Moustakides  and  S.  A.  Kassara 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Pennsylvania/ 

Systems  Engineering  Dept.' 
Philadelphia,  PA  19104 _ 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB,  Washington,  DC  20332 


14.  MONITORING  AGENCY  NAME  4  ADDRESS/'//  different  from  Controlling  Office) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Interim 


6.  PERFORMING  ORG.  REPORT  NUMBER 


6.  CONTRACT  OR  GRANT  N  UMBERf  sj 


AFOSR  77-3154 


to.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


61102F  2304/A6 


12.  REPORT  DATE 

1980 


1 3.  NUMBER  OF  PAGES 

Six 


t5.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 


t5«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 


DTIC 

ELECTE 
JUN  1  7  1980 


t7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  ft  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

Conference  on  Information  Sciences  and  Systems,  Princeton  University, 
March  1980;  paper  to  appear  in  proceedings 


t8.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Wiener  Filters,  Robust  Filters,  Spectral  Density 

tjj  T>h±  A’fe^r; 

~Z0.  ABSTRACT  (Contlnhe  on  reverse  side  It  necessary  and  Identify  by  block  number) 

Wu  Tonaidai1  the  Wiener  filtering  problem,  when  the  cros^  saectral_density 
matrix  of  the  signal  and  noise  is  not  exactly  knowifc  *  Vo  *o»taAn»f liters  whi 
are  saddle  point  solutions  for  the  criterion  of  pe^ 

MSE)  over  tji.^  classes  of  allowable  density  matric 
classes  are  given, 

*  /V 


which 

rm|uicefp(mean  square  error, 
Solutions  for  various 


D 


DD 


'T-n  1473 


EOITION  OF  I  NOV  «5  IS  OBSOLETE 


UNCLASSIFIED 


«irrnp,TV  r*  *4*»F|r»TlON  OF  THIS  PAGF  fWhen  Deta  Knteced) 


AEQSR] 


■% 


j  i  .i- —  — — -  * - — ■ 

NER  FILTERS  FOR  .CORRELATED  J 
IGNALS  AND  NOISEV  n - — 


.  r  a  / 


ABSTRAC  {Jts 

We  consider  the  Wiener  filtering  problem,  when 
the  cross  spectral  density  matrix  of  the  signal  and 
noise  Is  not  exactly  known.  We  obtain  filters 
which  are  saddle  point  solutions  for  the  criterion 
of  performance  (mean  square  error,  MSE)  over  the 
classes  of  allowable  density  matrices.  Solutions 
for  various  classes  are  given. 

I.  INTRODUCTION 

In  classical  Wiener  filtering,  we  need  to  as¬ 
sume  exact  knowledge  of  the  spectral  density  ma¬ 
trix  of  signal  and  noise.  In  many  applications 
this  assumption  of  exact  knowledge  Is  unreasonable. 
A  more  realistic  assumption  Is  that  our  matrix  be¬ 
longs  to  a  class  of  density  matrices.  This  class 
can  be  defined  according  to  our  knowledge  of  the 
true  density  matrix.  For  this  problem  we  will  de¬ 
rive  a  filter  that  performs  In  an  optimum  way  over 
the  whole  class. 

Nahl  and  Weiss  [1,2]  derived  the  bounding  fil¬ 
ter.  This  filter  Is  a  Wiener  filter  for  some  den¬ 
sity  matrix  Db  and  if  it  Is  applied  to  any  matrix 

from  the  class  the  MSE  error  is  bounded  by  the  min¬ 
imum  MSE  for  D^ .  Because  does  not  usually  be¬ 
long  to  the  class,  there  is  no  matrix  in  the  class 
that  can  reach  the  bound  of  the  MSE.  This  means 
that  there  is  a  possibility  for  better  performance. 

Kassam  and  Llm  [3]  derived  the  robust  Wiener 
filter,  when  signal  and  noise  are  uncorrelated 
(density  matrix  diagonal).  This  filter  sets  a 
bound  on  the  MSE  and  the  bound  can  be  reached  by 
some  matrix  from  the  class.  In  [4]  and  [5],  Poor 
generalized  some  of  this  work.  The  present  paper 
extends  the  above  idea  to  the  correlated  signal  and 
noise  case. 

Presence  of  correlation  is  possible  in  many 
applications.  An  example  Is  a  multipath  channel 
with  a  strong  signal  component,  weak  unwanted  mul¬ 
tipath  signals  and  regular  noise.  The  total 
"noise"  la  correlated  with  the  signal. 
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processes  are: 


i.  Dg (w) ,  Dn(v)  real,  even,  non-negative  functions 


^sn^ 


S  D  (w) » D  (w) 
s  n 


So  D  is  non-negative  definite,  and  diagonal  elements 
are  even  functions  of  w. 

Given  random  processes  s(t)  and  n(t)  with  den¬ 
sity  matrix  D  and  a  filter  h(t)  with  Fourier  trans¬ 
form  H(v) ,  the  MSE  for  signal  estimation  using  this 
filter  is 


e(D,H) 


:|s(t)-  / 


h(v)x(t-v)dv 


=  R  (0)-2  J  h(v)R  (v)dv  +  JJh(v)h(u)R  (v-u)dvdu. 


where  x(t)=s(t)+n(t)  and  R^^ (t )  is  the  cross  corre¬ 
lation  between  1  and  j.  Using  Fourier  transform  and 
Parseval's  theorem,  (2)  can  be  written  as: 

00 

e(D,H)-i- /  [Da(w)-2H(w)Dsx(w)+|H(w)l2Dx(w)]dw  (3) 


The  optimum  filter  for  D  is  given  by 


D8X(w) 


If  we  substitute  (4)  into  (3)  we  have  that  the  op¬ 
timum  MSE  is  given  by: 


VD)  ’  27 


Ds(w)Dx(w)- 


[Dsx(w)| 


II.  ROBUST  WIENER  FILTERS 

Lat  us  assume  that  our  processes  are  real, 
wide  sense  stationary  and  zero  mean  and  the  noise 
process  Is  additive.  We  also  assume  that  there  ex¬ 
ists  a  spectral  density  matrix  D  for  the  signal, 
noise  processes,  given  by: 


D(w)  D  (w) 

8  sn 

D*  (w)  D  (w) 
ill  a 


where  s  is  for  signal,  n  for  noise  and  (*)  for  com¬ 
plex  conjugate.  The  properties  that  characterize 

«Thls  research  Is  supported  by  the  Air  Force  Office 
U8AF,  under  Grant  AFOSR  77-3154. 


2.1  Definition  of  Robust  Filter 

Assume  that  a  class  &  of  density  matrices  Is 
given.  A  robust  filter  Hr  is  defined  by  the  fol¬ 
lowing  properties: 

a.  Hf  Is  an  optimum  filter  for  some  matrix  Drc&, 

so  that  e(Dr,Hf)  i.  e(Dr,H)  for  any  filter  H. 

b.  For  any  DcAwe  have:  e(D,Hr)  s  eop(Dr)*e(Dr,Hr) . 

Combining  a  and  b  we  have  the  saddle  point  relation 

•  e(D,Hr)  *  «op(Dr)  S  e(Dr ,H)  (6) 
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for  any  De  A  and  any  filter  H. 

2.2  Theorem  0. 

Let  A  be  a  convex  class  of  density  matrices. 
Then  the  pair  (Dr,Hr)  Is  a  saddle  point  solution 

for  MSE  over  class  A  and  the  class  of  all  linear 
filters  (It  satisfies  (6)),  If  and  only  If: 

V“r>  ■  “a  eoP(D) 

The  proof  Is  given  In  appendix  A. 

'  Theorem  0  Is  the  key  point  In  our  search  for 
the  robust  filter,  because  based  on  It  we  have  on¬ 
ly  to  maximize  e  (D) . 

op 

2.3  Maximization  of  e  (D) 

- op  — 

From  (5)  If  we  write  the  error  In  terms  of  s 
and  n  we  get: 


«ob(D) 

op 


2n 


Ds  *I>n<w)-|Dsn(w)|2 

D  (w)+D  (w)+2Re[D  (w) ] 

8  D  sn 


dv 


V) 


For  every  w  and  for  given  Dg(w) ,Dn(w) , |Dgn(w) | ,  the 
worst  RetDgn(w)]  is  "lDsn(w)l»  because  It  minimi¬ 
ses  the  denominator.  This  gives: 


S  (D) 
op 


»  (w)-D  (w)-|D  (w)j  2 
s  n  sn 

DB(w)+Dn(w)-2|D8n(w)| 


dw 


(8) 


For  given  Dg(w),Dn(w)  expression  (8)  as  a  function 

of  | 0  (w) |  is: 

sn 

a.  Increasing  for: 

0  £  |Dgn(w)|  t  mln{Dg(w),Dn(w)} 

b.  decreasing  for: 

■in{Dg(w).Dn(w)}  i  |DgJ1(w)|  <  /Dg(w),Dn(w) 

From  condition  a  and  b  we  can  say  Chat,  given 
Dg(w),  Dg(w),  the  worst  |Dgn(w)|  Is  the  one  that 

la  aa  close  as  possible  to  mln(Dg(w) ,Dq(w)} .  The 
worst  Dgn(w)  is  Dgn(w)  *  -lD,n(w)| • 

III.  APPLICATIONS 

1.  P,(w).  Pn(w)  given  with  upper  and  lower  bounds 

Let  0  £  L(w)  £  |D  (w)j  £  U(v),  with  L(w)  ,U(w) 
given.  If  we  define: 

A(w)  •  min{Dg(w) ,Dn(w)}  (9) 

then  the  worst-case  characteristic |d'(w) ! Is  given 

hr  “ 


_  f  L(w>  If  A(w)  s  L(w) 

|D._(w)j  «  J  A(w)  If  L(w)  i  A(w)  £  U(w) 

[  U(w)  If  U(w)  s  A(w) 

where  A(w)  defined  In  (9). 

From  section  2.2  we  liave  that  Dr  (w)“-| D  (w)l . 

sn  sn 

Figure  1  Illustrates  this  case. 

2 .  Upper  and  lower  bounds  on  (w)  and  Dp (w) ■ 

He  assume  that  bounds  L, (w) ,U, (w)  for  D, (w)  are 
given:  11  1 

0  £  Lj(w)  £  Dj^w)  sU^w)  l-s,n. 

In  addition  we  will  assume  knowledge  of  the  total 
power  of  signal  and  noise, 

m 

|oi(w)dw  «  2ito|  i»a,n 

m 

where  are  known. 

A.  If  there  are  no  bounds  on  |Dgn(w)| >  It  can 
reach  the  value  min{Dg(v) ,Dn(w)} .  Under  this  con¬ 
dition  the  error  Is: 


e  (D) 
op 


1_ 

2s 


mln{D  (w) ,D  (w)}dw 

s  n 


(10) 


Because  of  the  power  constraints,  there  are  several 
sub-cases.  First  we  will  give  some  definitions. 


Ag(w)  -  mln(Ug(w),  max[Lg(w),  Ln(w)]} 
Aj^(w)  -  min(Un(w),  max[La(v),  Ln(w)]} 
Bg(w)  -  min(Ug(w),  max[Lg(w),  Ug  (w)  ] } 


*n(w)  -  min{Un(w),  max[Ln(w) ,  Ug(w)]} 


Figure  2  Illustrates  the  definition  of  Ag(w) ,Bg(w) ■ 


Al. 


|Ag(w)dw  >  2so2,  (a 


,(w)dw  >  2so  ,  A  (w)dw  >  2»o 
n  n 


then: 


r  fLs(w)  if 

o'(w)  -  4 

*  [  l  (w)  otl 


Ag(w)  -  Lg(w) 


otherwise 


d!(w) 


Ln(w)  if  An(w)  -  Ln(w) 


t_(w)  otherwise 
n 


,_r.  2.2  1 
•op(D  )’°.+V2¥ 


(Lg(w)+Ln(w)-min(Ln(w) ,Lg(w)))dw 


Where  tg(v),  tg(w)  are  arbitrary  functions  with 


t,(w)  £  Ag(w)  and  la(w)  £  An(w) 


but  enough  for  Dg(v)  and  Dg(w)  to  fulfill  the  power 
constraints. 

"Oh.  .  o.iJ/Of 


m 


i 


□  o 


j!a(v)dw  >  2*o*,  |An(w>dw  <  2»a* 


D  (w)  is  as  in  case  AX 


1,0 

t  othe 


L(w)  -  A  (w) 
n  a 


otherwise 


where  i  (w)  arbitrary  function  with  t_(w)  i  A  (w). 
u  n  n 


.(O') 

op 

-  min(Lg(w), 

A3. 

m 

f  9 

m 

jBa<w)dw  >  2*o*  > 

jAs(w)dw 

* 

■yO* 

m 

m 

B  (w)dw  >  2 no2  > 
)  a  n 

|An(w)dw 

-40 

also  assune 


^W°n  -  fAn(w)  >  2*a*  -  Jab(w)(1w 


’  Aa(w)  if  As(w)  -  Ba 

Ds(v)  ”  *8(v)  otherwise 

A  (w)  if  A  (w)  -  B  (w) 
D>)  “  n 

t  (»)  otherwise 

* 


^  otherwise 

Where  Ag(w)  S  la(w)  <  Bs(w)  and  Dr(w)  £  ift(w) 

S  *n(w) 

«op(Dr)  •  a*  -  ^|{La(w)-«in(Un(w),La(w))>dw 

m  ~  m  m 

AA.  |Ba(w)dw  >  2*o*  >  |Aa(w)dw,  2*o*  >  |Bn(w)dw 


D‘(w)  as  in  case  A3. 

fu  fa)  1* 

D'(w)  -  4  B 

[t  M  ot! 


0  (w)  -  A  (w) 
n  n 


otherwise 


ia(w)  >  D*(w) 


•^O)*)  -  o*  -  ^  pLa (w) -min  (Un  (w)  ,La  (w)  )  }  dw 
2*0*  >  |«a(w)dw,  2»o*  »  j^fajdv 


Dr(w)  J  °«>(W)  “  "»<W)  ■  B- 

*  ]  »a  otherwise 

f  V  (w)  if  0  (w)  -  B (w) 

r  I 

|  in(w)  otherwise 


otherwise 


Where  i  (w)  >  B  (w)  and  l  (w)  >  B  (w). 
8  8  an 


•  ST  jj 


nln{lla(w),  Un(w))dw 


In  the  cases  above,  interchanging  the  roles  of  s  and 
n  in  the  conditions  give3  us  similar  results.  The 
proof  for  cases  A1  and  A3  can  be  found  in  Appendix 
B. 

B.  Suppose  we  are  given  a  function  R(w)  such  that 
0  a  |Dat|(w)|  £  »(w)  <  min(La(w),  Ln(w)) 

Because  of  the  maximization  problem  |Dsn(w)|has  to  he 
as  close  as  possible  to  mln(Da(w),  Dn(w)>,  so 

I  Dan(w)|  •  R(w).  To  find  the  pair  Da(w),  D*(w) 
notice  that: 


e  CD) 
op 


x  Ds(w)-Dn(w)-R  (w) 

2*  D  (w)+D  (w) -2R(wT 

—  8  n 


Ob 

JL  +  i  S(w)  N(w) 

2*  R(w)d  +  2*  S(w)+N(w) 


Where  S (w)-Dg(w)-R(w)>0  and  N (w)-Dn(w)-R(w)>0.  To 

maximize  (11)  it  is  enough  to  maximize  the  second 
berm.  But  this  term  is  the  expression  for  the  min¬ 
imum  MSE  for  uncorrelated  signal  S(w)  and  noise 
N(w),  with 

La(w)  -  R(w)  i  S  (w)  £  l>a(w)  -  R(w) 

L  (w)  -  R(w)  £  N(w)  <  0  (w)  -  R(w) 
n  n 


•ad  |S(w)dw  •  2na 


]H(w)dw  -  2iro 


jR(w)dv 

m 

[R(w)dw. 


This  problem  has  been  solved  in  [3]  and  gives  the 
solution  to  the  present  case. 

3.  Given  classes  for  D^(w) ■  D(w) 

Ws  assume  again  knowledge  of  the  total  power 
of  signal  and  noise.  When  there  is  no  restriction 

Oft  D _ (w)  we  have  seen  that  the  error  is  given  by 

•a 

(10)  and  It  can  be  written  in  the  following  way: 

r 


i 


vhe^e 


g(x) 


[  x  for  0  <_  x  <_  l 
1  for  x  >  1 


Da(w) 


-*  |2 


1 2 


•-«x-  l*J  '“mI 

— <  (1-e)  -22-  +  c  -22- 
D*  —  D'  D" 

X  X  X 


i2 


D" 


|(1-€)d;x  +  eD»J2  <  ^(l-c)l  |t  +e — g5  }  . 


g(x)  is  a  convex  function  of  x  for  x  a  0.  So  (12) 
can  be  used  as  a  measure  of  distance  between 
D#(w)  and  Dn(w).  Now  for  given  classes  of  Dg(w) 

and  &n(v)  we  want  the  pair  Dg(w),  D^'(w)  that  has 

the  suxlmuo  distance.  It  turns  out  as  we  can  see 
from  [4,6]  and  [7]  that  this  pair  does  not  depend 
on  the  form  of  g(x).  As  long  as  g(x)  Is  any  con¬ 
vex  function  the  pair  Is  always  the  same.  Solu¬ 
tions  to  this  problem  are  also  given  In  these  works 
for  various  classes. 

IV.  NUMERICAL  RESULTS 

He  will  assume  the  case  when  signal  and  noise 
are  given  as  in  Figure  3  and  there  is  no  restrict¬ 
ion  on  |D  (w) | .  The  robust  filter  is  simply: 


Bjw)  -4 
.  1° 


when  D  (w)  >  D  (w) 
s  n 


0  otherwise 


It  turns  out  that  this  filter  has  the  same  error 
performance  for  any  D^w).  If  we  use  instead  of 

Hr(w)  the  filter 


Hu(w) 


Dg(w) 


Ds(w)+Dn<w) 


assuming  that  signal  and  noise  are  uncorrelated, 
then  the  error  is  given  by  (3)  and  It  becomes  max¬ 
imum  for  _  ,  . 

D  (v) 


an 


->'Ds(w)Dn(w). 


0  i.  c  s  1  then  the  expression: 


0<*,w) 


D*(v)-D*(w)-|0^x(w)r 

0‘(w) 


is  a  convex  function  of  c. 

Proof.  It  is  sufficient  to  prove  that: 

C(c,w)  a  (l-c)C(0,w)  +  tC(l,w). 

Subtracting  each  side  from  D*(w)  we  have  to  prove 
that 


♦{ (l-e)D^  +  eD"  ) 


But  this  Is  the  Schwartz  Inequality. 

Proof  of  Theorem  0. 

The  only  If  part  is  easy.  From  (6) 

e(D,Hr)  £  e(Dr,Hr)  -  eop(Dr) 

but  e(D,H  )  i.  e  (D)  so  e  (Dr)  i  e  (D) . 
r  op  op  op 

To  prove  the  if  part,  define  De  •  (l-e)Dr  +  eD, 

0  £  t  <  1  where  De  A.  Because  of  the  lemma  and  the 
fact  that 


eop»E>  ’ 


w)dw 


e  (D  )  is  a  convex  function  of  e.  So: 
op 

e  (De)  >  (l-e)e  (Dr)  + e  e  (D)  and 
op  op  op 


o  x  -St 


e  (De)  -  e  (Dr) 


op 


i  e  (D)  -  e  (D  ) 
op  Op 


(13) 


Because  -2E. 


e _ (DE)  -  e  (Dr) 


is  monotonic  with  respect 


In  Table  1  are  given  some  numbers  for  differ¬ 
ent  b.  The  error  is  er  when  we  use  Hr(w)  and 

la  the  worst  error  when  we  use  Hu (w) .  In  the  third 
column  we  can  see  the  percentage  of  performance 
Improvement.  Also  shown  is  the  optimum  error  e*  of 
H  <w)  when  D  (v)-0. 

APPENDIX  A 

Before  proving  Theorem  0,  we  will  prove  the 
following  lemma. 

1.  Let  D\  D"cA,  define  De-(l-t)D'+  eD" 


to  e  and  bounded  as  we  can  see  f |om  (13)  its  limit 
exists  as  e-»-o+.  But: 

Vp(D  )~eop(D  )  .  1_  |  dw  (U) 
Now  ^  G(c,w)  is  also  monotonic  with  respect  to  r. 
From'  (13)  and  (14) ,  we  have 


0  i 


11m 

e*o 


mil  G(e.w) 

2* 


dw 


1  11m  G(e ,w)~  G(o.w)  . 
2r  e-*-o+  e 


m 

G(e,w)lE.0+ 


dw 


But  the  last  expression  after  we  Cake  the  deriva¬ 
tive  and  make  some  manipulations,  equals  e(D,Hr)  - 

(Dr).  So  e(Dr)  t  e(D,Ht). 


op' 
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For  the  second  application  we  will  outline  the 
proofs  for  subcases. A1  and  A3.  In  a  similar  way, 
*e  can  prove  the  rest  of  the  cases. 

Umm  2.  If  a,b,c  non-negative  numbers  with  a  >  c 

then: 

a  -  c  i  mln(a,b)  -  mln(fc,c) 

Proof.  a  -  c  ^  0  i  mln(a.b)  -  k 

,  a  -  c  »  mln(a.b)  -  c  so 

a  -  c  i.  min(a.b)  -  min(4,c) 

From  (10)  we  have  that  the  error  is 

m 

e  (D)  -  -yjimin{D  (w)  ,D  (w)}dw 
op  2iu  s  n 

Case  Al. 

*ln{D  -D  )  -  mln(L  ,D  }  + 
s  n  s  n 

(min (D  ,D)  -  mln(L  ,D  }] 

•  n  8  n 

(using  Lena  2)  i.  min{Ls,Dn)  +  Ds  -  Lg 

(using  Lemma  2)  s.  min{L  .L  }  +  D  +D  -L  -  L_ 

8  n  s  n  8  n 

So 

•op®5  *  °s  +  °n  "  ^flVw,+Ln(w)- 
mln{Lg(w) ,Ln(w)} ]dw 

And  we  have  equality  when  D  ■  Dr. 

Case  A3. 

mint D  , D  )  *  mln{L  ,D  }  +  [mln{D  ,D  }  - 
8  n  s  n  s  n 

mln{Ls,Dn)]  s  min{Lg,Dn>  +  Dg  -  Lg 

i  mln(L  ,B  )  +  D  -  L  and 

s  n  s  s 

m 

VD)50^1  [  Lg (w) -min( Lg (w) ,Un(w)} ]dw 
with  equality  when  D  *  Dr. 
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TABLE  1 


b 

e 

r 

eu 

Z 

ei 

0.10 

0.09 

0.12 

34.5 

0.07 

1.00 

0.50 

0.62 

25.7 

0.33 

5.00 

0.83 

1.10 

31.7 

0.62 

10.0 

0.91 

1.20 

34.5 

0.73 

TABLE  I  Performance  Comparlslon  of  Robust  and 

Nominally  Optimum  Filter  for  Example  of 
Figure  3.  (No  bound  on  |Dr  (w) |  ). 


